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Abstract 

We employ the Riemann-Hilbert problem for solution of the initial-boundary value problems for nearly 
integrable equations on the half line which have important applications in physics. The detailed derivation of 
the integrable and perturbation-induced evolutions of the spectral data are given for the two most important 
models, the self-induced transparency and the stimulated Raman scattering. In particular, we prove that 
dephasing stabilizes the spike phenomenon in the stimulated Raman scattering. 

1 Introduction 

Nonlinear integrable models of resonant interaction of radiation with matter possess a number of features which 
distinguish them from a variety of integrable systems. First of all, these equations determine initial-boundary- 
value problems, because, along with an initial condition, some boundary conditions are necessary to properly 
formulate the problem. Secondly, they have singular dispersion relations [jl], As examples of such integrable 
equations we mention the Maxwell-Bloch equations of self- induced transparency equations of stimulated 
Raman scattering ||] , the Karpman-Kaup equations 0, ^) , the equations of the interaction of the electrostatic 
high-frequency wave with the ion-acoustic wave Q, and the equations of the resonant interaction between 
optical solitons and impurity atoms in fibres || . 

There exists a well-developed formalism for solution of the integrable equations with the singular dispersion 
relations when the boundary conditions are asymptotic, i. e. they are posed at the plus or/and minus infinity 
[H [u| O. Along with the fact that the semi-infinite interval, when one of the boundaries is placed at a finite 
point, e. g. zero, is physically more appropriate, the formalism of the inverse scattering has to be substantially 
modified as compared with the whole line case |0, O . 

Recently significant progress has been achieved in solving the nonlinear evolution equations with singular 
dispersion relations posed on the semi-infinite |l4], [l5| and finite |16, [lTj intervals. Also it was realised that 



the Riemann-Hilbert (RH) problem is a natural setting for analysis of the above nonlinear models. The data 
uniquely characterising the solution to the RH problem serve as parameters of the solution to these nonlinear 
equations. 

All the above-mentioned papers deal with the integrable evolutions. On the other hand, in the experimental 
observations one has various additional (frequently small enough) effects which destroy the integrability property 
but must be accounted for complete explanation of the experimental results. For instance, the phenomenon of 
the spike formation in the pump depiction zone in the stimulated Raman scattering |lS[ ] is directly related with 
the dephasing effect in the medium [|l9|. Hence, a problem arises to describe evolution of the RH (or spectral) 
data in the presence of small perturbations, i. e., to develop a perturbation theory on the half line. 
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In the present paper we extend the perturbation theory developed for the whole line [|20| |2l|], to the half line. 
We take the equations of the self-induced transparency (SIT) and stimulated Raman scattering (SRS) posed 
on the half line as the examples. In both cases the dephasing effect, which destroys the integrability of the 
models, is considered as the perturbation. We show that for the half line, as for the whole line, construction of 
the perturbation theory amounts to introduction of an additional matrix functional, which we call the evolution 
functional. Moreover, it turns out that the evolution functional has the explicit expression very much similar 
to that in the case of the whole line (see, for instance 21 ). 

The plan of the paper is the following. In section 2 we briefly remind the formulation of the initial-boundary 
value problems for the SIT and SRS equations. In sections 3 and 4 we present necessary results on the direct 
and, respectively, inverse spectral problems for the half line. In particular, we demostrate the advantage of 
dealing with the analytic functions with constant determinant for derivation of the evolution equations for 
the spectral data in the case of the half line. Our main section is section 5, where the perturbation-induced 
evolution equations for the RH data are derived. 

This paper is dedicated to Professor Pierre Sabatier whose versatile activities have determined in many 
respects the contemporary status of inverse methods in science and technology, in honor of his 65th birthday. 



2 Nonlinear models with singular dispersion relations 

The SIT phenomenon consists in a soliton-like 27r-pulse propagation in a medium without energy loss and is 
described by the Maxwell-Bloch equations Q 

£ x =<p>, Pt + 2iXp + Tp = N£, N t = ~(£p + £p). (1) 

Here £(x,t) is the envelope of the light-pulse electric field with the carrier frequency close to the particular 
transition frequency of an atom, 

"-(? £ 

is the density matrix of the atomic ensemble (medium) with p(x, t, A) being induced polarisability and N(x, t, A) 
being normalised population difference between the chosen atomic levels. We have detM = — 1 and for atoms 
initially on the lower level N(t = 0) = — 1 and M[t = 0) = 03. A is the normalised frequency mismatch between 
the carrier frequency and atomic transition frequency, T is the dephasing factor. We think of the value of T 
as being small, r <C 1, in order to justify the applicability of the perturbation theory. The coordinate x is 
the resonant medium extension and t is time. The symbol < . . . > means the averaging over the frequency 
mismatch with the normalised distribution function <?(A): 



<p>= dXp(x,t,X)g(X), / d\g(X) = 1. 
J —00 J —00 

As regards the SRS equations, we follow the approach by Leon and Mikhailov JlJ]: 

u t + Tu=-g dkax^e^ 21 ^ , a lx = ua 2 e 2lkx , a 2x — -uaie~ 2ikx . (2) 



Here ai(x,t, k) and a,2(x,t,k) are the electric field envelopes of the pump and Stokes pulses, respectively, the 
parameter k describes a spectral extension of the input pulses, u(x,t) is related with the polarisability of a 



2 



Raman-active medium, g is a coupling constant. It should be noted that the parameter fc can be connected 
with the group- velocity dispersion of the pulses in the medium H] . 

Both nonlinear systems are considered in the quadrant 0<a;<oo,0<i<oo with the initial and boundary 
values: 

a) SIT: £(0,t) = E (t), M(x, 0, A) = M = ( ~ N ° ) , £^0 as t oo; (3) 



Po iV< 



b)SRS: u(ar,0)=0, ai(0,i,fc) = h(t,k), a 2 (0, t, fc) = J 2 (t, fc), (4) 

Note that for the SIT equations we have an initial-boundary-value problem, while for the SRS equations we 
consider a pure boundary-value problem. 

For r = both systems (1) and (2) are integrable and admit the Lax representation ||, ^2|. In both cases 
the dispersion relation is non-analytic (singular) function on the complex plane of the spectral parameter. 



3 Direct spectral problem 

Because both the SIT and SRS equations are solved by means of one and the same Zakharov-Shabat spectral 
problem j2^], we will analyze below the direct spectral problem for both systems together (omitting below the 
explicit dependence on the other coordinate). The spectral problem is written as 

X£ = -ik[cr 3 ,x] + Qx- (5) 

Here x is t ne 2x2 matrix function, 

q 
-q 

is a potential with £ = 2q for the SIT and u — q for the SRS. The coordinate £ stands for t in the case of the 
SIT and for x in the SRS case, fc is a spectral parameter coinciding with the mismatch parameter for SRS. It 
should be stressed that the physical status of the potential q is different for the SIT and SRS. Whereas for the 
SIT the potential q (or £) is the main experimentally measured quantity, for the SRS the potential u could be 
reconstructed only indirectly from the measured intensities |ai| 2 and |o,2 1 2 - 

Let us define the matrix Jost solutions J±(£, fc) of the spectral equation (5) by means of the following 
conditions: 

J_(0,fc)=7, J + (£,k)^I, £^co, (6) 
where I is the identity matrix. The scattering matrix is defined in the usual way: 

S = E- 1 J+ 1 J_E, E{£,k) = e-^ 3 , lmfc = (7) 

and has the structure 

a —b 



S=\l _ ) , ao + 66=l. (8) 

The scattering matrix is connected with the boundary value of J+, S(k) = J^I 1 (0, fc). Note that the potential 
is anti-Hermitian: Q* = —Q. Hence, we have the involutions (Imfc = 0): 

4&*) = J±\t:k), S\k) = S-\k). (9) 
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The spectral equation (|5|) can be represented in the integral form: 



j(£,fc) = £(0 



j + /' d?Er i (OQ(?)j(e,k)Em 



E~HZ), (io) 



where Jo = J(£ = Co)- The Jost solution J_ ( J + ) corresponds to the choice £o = (£o = oo). It follows from ( |To| ) 
that the matrix functions (J_ ,J+ ) and ( J^ 1 , J_ ) are solutions to equation (^J) and moreover holomorphic 
in k in the upper and lower half-planes, ±Im/c > 0, respectively. Here J^ n ' stands for the corresponding column 
vector. By noticing that a solution to (0) transforms to another solution if multiplied on the right by a diagonal 

^-independent matrix and that det(J_ , J+' ) ) = a and det(j| 1 ' ) , ji 2 ') = a wc introduce the solutions to (||) 
having unit determinant: 

^(J^a-Vf), *_ = (a" 1 jf J ji 2 )). (11) 

By definition the functions k) and ^-(^, k) are analytic in the upper and lower half-planes of the complex 

fc-plane, respectively. On the real line Imk — they can be given in terms of the matrix Jost solutions and 
entries of the scattering matrix: 

= J+EG+E- 1 = J_EG-E~ l , = J+EH^E- 1 = J_EHZ 1 E- 1 . (12) 

Here 

G ( a \ G __ ( 1 b/a 



b 1/a J ' \ 1 

H+ = ( l/o ) ' ^ = ( 6/5 1 
As follows from (^|) and (|l2|), the matrix functions ^±(S,,k) satisfy the involution: 

tf + (k) = ^Z 1 (k). (13) 

The representation ( |l2] ) involves the factorizations of the scattering matrix with division: S = G+GI 1 and 
S = HZ 1 !!-. The advantage of using such factorizations for the half line stems from the fact that both G_ and 
H— contain a single quantity: b/a and b/a, respectively. On the other hand they provide the values of ^± at 
the boundary: f+(0, k) = G_ and *_(0, k) = HZ 1 . In view of this, it is convenient to introduce the following 
notations: 

m = ^, ai > k ) = 4ry ( 14 ) 

a(k) a(k) 

The functions (3(k) and a(k) are usually referred to as the reflection and transmission coefficients, respectively. 
Then the boundary values of the solutions v I / ±(^, k) read 

mc=o)=(; i), Me=o)=(_^ as) 

For completeness, let us present also the asymptotic values of \&±(£, k) as £ — > oo: 



l/a \ / a -/3/ae- 2lfc « 



a ), o " 1/fi ) (i6) 
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Now, let us consider in more detail the analyticity properties and the fc-asymptotics of the matrix functions 
^i. These properties are easy to deduce from the integral equations for these functions: 



*+i 2 \ = /0\_ / jf dCe- 2 ^-^ 



^22 



*-i2 \ _ f \ / / ? dC'e- 2lfc («-«')gvl/ 



-22 



*_ 22 ~ 1 ^ _ f€ 



Evidently the first column Vl/V/ and the second column vff^ are entire functions of k, while and 
are meromorphic functions in the upper and lower half-planes of the complex fc-plane, respectively (see also 
Ref. By application of ( p"5| ) we conclude that the reflection coefficient /3(fc) is meromorphic in the upper 

half of the complex plane. Moreover, as 0/a = b, by using the identity J + i2(0, k) = b{k) (see (0) and (||)) we 
conclude that the quotient (3{k)/a(k) is holomorphic in the upper half plane. Equivalently, a(k) and b{k) are 
holomorphic in the upper half of the comlex fc-plane. Note that the poles of /3(fc) coincide with zeros of a(k). 

From the above integral equations it is seen that ^>± — * / as k — ► oo. Integration by parts in the integrals 
containing exp(±2ifc£') gives the generalized asymptotic expansions: 

*+2i J \0 J 2ik\ g(0 - e 2lfe «g(0) / v> 2 

?(0 A , o ( l 



( * + 12 X 




! 




V *+22 , 


H 




+ 2ik [ 



\/0\ J_/ ?(0 - e- 2 ^ g(0 ) \ /l 

*-22 J U / 2i* V -/o d^'k(OI 2 + <z(0)J * dC'e- 2lfe «'g(^) J + U 2 

These formulae give triangular matrices for ^±(0, k) in agreement with the the boundary values (|l5|). For 
Imk = and £ 7^ 4'±(£, fc) do not have asymptotic expansions in the usual sense due to the presence of the 
terms exp(2i/c£)g(0) and exp(— 2ik£)q(Q). 

Now we are in a position to express the potential Q through the matrix elements of SSf±: 

g(0 = 2i lim fc*+i 2 (£, fc), = 2i lim fctf_ 2 i(£, *). (17) 

k— >oo /c— >oo 

Alternatively, one can use the following formulae: 

q(0 = 2i lim k (*_i2(£ fc) + e- 2lfc «/3(A:)) , 
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q{£) = 2i lim fc (* +2 i(fc k) - e 2ik ^(k)) 



Here we have taken into account ( |l7| ) for £ = with , 3> + i2(0, fc) = (3(k) and 4'_2i(0, fc) = — /3(fc) from (|i"5|). 

The reconstruction of the potential requires knowledge of the functions ^>±. The latter can be obtained via 
solution of the RH problem, which we consider in the next section. 



4 The Riemann-Hilbert problem 



We have from © 
where 



y + = y-EGE-\ lmfc = 
G(k) = H_G_ = H+G+ -- 



1 (3 
P 1 + I/3I 2 



(18) 
(19) 



This is the matrix RH problem, i. e., the problem of analytic factorization a nondegenerate matrix G(fc), given 
on the real line Imfc = 0, into a product of two matrices analytic in the upper and lower half-planes, respectively. 
The above RH problem has the canonical normalization condition: 



I 



as k — > oo. 



(20) 



To completely characterize the solution to the RH problem (|18[)-([20[) one needs to specify the spectral data. 
First, it is the value of j3(k) on the real line Imfc = 0, which enters G(k) Jl9| ) and represents the continuous 
part of the spectral data. Second, the discrete spectral data should be given. The latter include zeros, kj, 
Imkj > 0, of the function /3 _1 (fc) [(3~ 1 (kj) — 0] and the transition coefficients Cj. The transition coefficients are 
given by the residues of ^ + or at the points kj and kj, respectively. Using the identity det v E'±(^, fc) = 1, 



the involution (13) and formula (p_lj) it is not difficult to obtain (see also J14|) 



Res 



*f(fc), 



kj =Cj^ 



(i) 



(kj), Res 



^\k),hj = 



-Cj® 



(2) 



(kj)- 



(21) 



The ^-dependence of c,- follows from the equation 



0« 



-ik (cr 3 + 1) 



Res^\k 3 



+ Q 



Substituting here @, we get Cj(£) = c} 0) e" 2lfc ^. 

The last step towards the complete description of the RH data involves the analysis of their temporal 
behavior. Here we analyse in detail the two models of section 2, the SIT and SRS equations. These models are 
integrable and hence represent the compatibility condition of a linear system, i. e., the Lax pair. The first Lax 
equation is the spectral equation (||). The time evolution is determined by the second Lax equation: 



a r *± = y±En ± E- 1 + V±V±. 



(22) 



Here t — x for the SIT and r = t for the SRS models. The matrices V±(£, r, k) are the limitimg values, as 
k — > Refc ± iO, of a piece-wise holomorphic function V(£, t, fc). The choice of V(£, r, fc) depends on the model. 
In the same way, the dispersion relation f2(fc, r) is also a piece- wise holomorphic fuction with the limiting values 
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Q±(k,r) defined for real k. (Hence the name "singular", i. e., non-analytic dispersion relation). The matrices 
V and satisfy the involutions 

Vl(k) = -V-(k), n1(k) = -Q-(k), 



which follow from those of ^(k) (|13|). Note that the compatibility of the Lax pair (||) and ( J22] ) is guaranteed if 

V is chosen in such a way that the jump matrix AV (£, r, fc) = V+ (£, t, k) — V- (£, r, k) satisfies 

d s AV(£, r, k) = [-ika 3 + Q, AV(£, t, k)], (23) 

while f2(fc, r) is still arbitrary. It means that the dispersion relation is not defined by the model. In fact, as it 
was realized in Ref. it is specified by the imposed boundary conditions (see also |l4|| ). Let us first obtain 

V for the models we consider. It is convenient to represent a V satisfying (|23j) in the form 

1 f°° d\ 

V{t,r,k) = — ——J-(Z,T,\)E(Z,\)AV(0,T,X)E- 1 (Z,\)JI 1 (Z,T,\), (24) 

11:1 J-oa A ~ & 

Formula ( pi|) trivially transforms in the SIT case by letting 

AV(£ = 0) = -^M , (25) 



while for the SRS equations we have 

"Lt ) = hJ W + J 2 J<V* = 7 X *W + I 2 *1V 4 * (26) 



thereby 



a 2 e 



Note that in both cases trV — trAV = 0. 

To derive evolution equations for the spectral data and the dispersion relation f2 we recall that the analyticity 
properties of ^±(k) depend on the conditions @ defining the Jost solutions to the spectral equation (||). 
Hence, to have analyticity of ^±(k) for arbitrary r one must demand that those conditions hold for all r. 
Equivalently, we can deal with the boundary and asymptotic values of X P±(^, t, k), which must satisfy ( |l5| ) 
and ( |l6| ) for arbitrary r. Substitution of formulae ( |l5| ) and ( |l6| ) into the second Lax equation (^2|) supplies 
one with the evolution of the spectral data and a together with the dispersion relation O. After somewhat 
lengthy but straightforward calculations (in computing the limits of V± as £ — ► oo one uses the identity: 
V.p.{exp(±2ik£)/k} -> ±i7T<5(fc)) we obtain (see also Ref. Q) 

/?r = -^x/? 2 + 2^/3 + 4i2 > (28) 
a T = ((SA^Wfif-iJaa _ ^ _ v io)^ a (2Q) 

Here V^j = (£ = 0) . For completeness we give also the dispersion relation: 



fi+ = 



py+21 v +u 



v +21 > JV +21 v +22 
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0V® 



y(0) 



v -12 



V 



(o) 



Here we note that the corresponding evolution equations for the entries (a and 6) of the scattering matrix 
are highly nonlinear coupled equations; while for the quotient (3 = b/a one gets a simple Riccati equation. This 
is the advantage of using the matrices i f?± defined by ([[l])- The relevance of the Riccati equation for equations 
with the singular dispersion relations was noted previously (see, for example, p3 for SIT and |l4|, |2(| for the 
SRS). 

The time evolution of the discrete RH data is found as follows. First, let us find the time dependence of Cj. 
It follows from (pl|) taken at £ = that 



thereby 



cf\r)=Re S [f3(k,T),k j }, 



c j (Z,T) = e- 2ik ^Bes[0(k,T),k j }. 



Evolution of the poles of /3 is derived in the following way. Inasmuch as (3 1 {kj 1 r) — for all r, we have 

dki {d/dr)l3- x {k) 



^3 

dr 



which, in view of ( p8| ) and the identity (d/dk)f3~ 1 (kj 



(30) 



(31) 



1/c) , gives 



dr 



= -^° 2 ) 1 (fc i )cf. 



(32) 



We remind that 



a) SIT: Vl°^[x,k 



b) SRS 



i 

' 1 
2i 



dA 



-<?(A)p (0) (x,A) 



dA 



h(t,X)h(t,X). 



Hence, kj = const for the "causal" solutions (in Zakharov's terminology J27j) for the SIT or in the absence of 
the Stokes pulse for the SRS. 

In this paper we are interested primarily in description of the effect of a perturbation on the evolution of the 
spectral data rather than characterization of specific solutions to the integrable nonlinear equations. Therefore, 
we proceed to the derivation of the perturbation-induced evolution of the spectral data. This is done in the 
next section. 



5 Perturbation-induced evolution of RH data 

Let us attribute the "variational" derivatives to the perturation-induced evolution. For instance, we introduce 
the perturbation matrix R as the variation of the potential, 
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A variation of Q leads to the corresponding variations SJ± of the Jost solutions. We have 

(J± l 5J±)t = -ik [a s , J± l 8J±] + Jg 1 SQJ±. 

Solving this equation we get 



SJ ± = J±E(0 ^ d^'E- 1 (e)J± 1 SQJ±E(C 



) E- 1 ^) 



with £o=0 for (5J_ and £0 = 00 for 8J + . Recalling the definition of the scattering matrix (|?]) and using the 
relations (|l^), we obtain its variation: 

A 

— = G+T+(fc)Gl 1 = H^T^H-, (34) 
or 

where T±(fc) = T±(0, 00; fc) and 

rp 

T±(g,p;fc)=/ d£E- 1 V± 1 m±E. 

J Q 



By variation of the relations (|12|) and using the above formulae one easily derives the perturbation-induced 
evolution of the matrices i B±: 

= V ± EIl ± E-\ (35) 

OT 

where the functional II is defined as follows (omitting the explicit r-dependence) 

TT (P V\- ( T +ii(°'^ fe ) -T+i 2 (C,oo;fc) \ u(fk ,_( T_ u (0,e;fc) T_ 12 (0,£;fc) \ 

- ^ T+2l(0 , efc ) -T +11 (0,£;fe) ) ' U "«'* J " ^ -T_ 21 (£,oo;fc) -T_ u (0,C;fc) J " W 

The r.h.s. of ( ^H|) should be added to the second (evolutional) Lax equation ( Bgh to account for the perturbation- 
induced evolution of \&±. Note that trll± = in agreement with the identity det\l/± = 1 and n^_(fc) = — II_(fe) 
due to the involution ([[si). 

Formula ( ]35| ) accounts for the perturbation exactly. However direct application of (|3^) is not possible due to 
the fact that *ff± should be obtained from the RH problem, solution of which requires knowledge of the spectral 
data, while the perturbation- induced evolution of the latter explicitly depends on \I/± (see below). On the other 
hand, if the perturbation is small, one can expand the exact evolution equations for the spectral data into series 
of approximate equations which are solvable. 



5.1 Evolution of the continuous datum 

In view of the boundary values of if? Jl5|), the perturbation-induced evolution (variation) of (3(k) obtains from 
( j36| ) taken at the boundary £ = 0. We have 

8/3/Sr \ ( 1 P\(0 -T+i 2 (fc) 
j \ 1 \ 



Hence the variation of (3{k) satisfies 

1*00 

2 (k) = - d£e 2ik t(^ 1 m + ) 12 . (37) 
Jo 



5(3 _ „ 
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Substitution of the expression for R from (p3|) gives 



5/3 



d£e 



2ifc£ 



r(* + i2) 2 +r(* +22 ) 2 



(38) 



The integrand in ( |38| ) is meromorphic for Imfc > by the definition of *$>±{k) (|ll]) and this is in agreement of 
the analytic properties of (3{k). 

Equation ( 38 ) allows us to prove the applicability of the linear limit in accounting for the effect of non-zero 
dephasing on the spike formation in the SRS equations pij . Moreover, the argument we present below is valid 
for non-small dephasing parameter T. Following the arguments of Ref. |24| (see also Ref. |l4]]) we note that 
the spike of pump radiation, first observed in Ref. (T^| (see also more recent Refs. [^8| and J29|), occurs at a 
time to = to when /3(fco,ro) = 0. Hence, the limit (3 <C 1 serves as an approximation when dealing with the 
spike formation. For small (3, solving the RH problem ( |l8| ) in the linear limit and using (|l7| ) we get (in the SRS 
notations (||)) 

u = -- [ dke" 2lkx (3{k) + 0{J3 2 ). (39) 



Here the integration is in the complex plane above all poles of (3(k). On the other hand, substitution of 
8u/5t — —Yu in equation ( |38| ) gives 



8J3_ 
~5t 



dxe 2lfe ^ + 0(/3 3 ) 
ot 



-r / dxe 2ikx u + 0(/3 3 ) 
'o 



-r/3 + 0((3 3 



(40) 



Here we have used formula (|39|). 

Equation ( ^0| ) is a quite remarkable result. Indeed, the integrable part of the evolution of /3(fc,r) is given 
by the Riccati equation (|28|), which is quadratic in (3. When (3 is small (spike formation), addition of the 
perturbation-induced evolution in the same order in (3 to the integrable limit requires keeping the terms linear 
and quadratic in (3. Howerver, in equation ( |40| ) the quadratic term is equal to zero. Thereby, the only effect of 
the dephasing term is stabilization of the spike of pump radiation in the SRS as it was realized, though without 
a rigorous argument, in Refs. [Q. 

5.2 Evolution of the discrete data 



Here we describe the action of a perturbation on the evolution of the discrete RH data. A perturbation-induced 



evolution of cf^ is calculated from (|30|) and (38). Namely, 



Res <i J do 



2ik(, 



f (* +12 ) 2 + r(* +22 ) 2 



, kj 



(41) 



Separating out the singular parts of the functions inside the integrand by using formula ([21 



#+12 (fc) 



#+ 22 (fc) 



c (0) e -2ife^ 



c (0) e -2ifc J « 



J 



*+n(fci) + *+f 2 (*0, 



y +21 (k 3 ) + ^l 2 (k), 



(42) 
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we can write equation ( pl| ) as follows 



5c 



(0) 



St 



= Res 



.(0)2 



(k-k])J 



dfe 



c (0) -oo 



(43) 



Here we have introduced a independent, and a regular function of k, -F 2 (£, k), by (omitting the explicit 

r-dependence for simplicity) 

Fi(0 = m (*+nK, fc,)) 2 + r(0 (* +a iK, ^)) 2 , 



Calculating the residues in (E3) , we obtain the perturbation-induced evolution of the coefficient c 



(0). 



6c 



(0) 

-1— - 2c 
5r J 



(0) 



/>OC />oo 

f d£F 2 (%,£) + icf / d^e- 2 *^F!(0 

JO 



(44) 



Finally, the perturbation-induced evolution of fcj is found as follows. We start with the formula similar to 
( fjl| ) but with the variational derivatives, i. e., 



5k 



j_ 

St 



{5/5T)(3-\k) 



Using the relations 



0-\k) = 



<D[(k - k 



{d/dk)0- x {k) 
dp- 1 



(45) 



■3) J> 



dk 



(k) = — + 0(k-k J ), k^k 



3' 



which follow from (BOf), and taking the perturbation- induced evolution of (3 in the form of (37), we rewrite (}45|) 



as 



5kj_ 
5t 



(fc-A: J ) 2 T +12 (fc)(l + 0(fc-fc J )) 



0(k - kj 



(46) 



Note that the numerator is non-zero because T + i 2 (fc) is meromorphic function having poles of the second order, 
as it is seen from (|42|). Using ( |38| ) and ji^), we obtain 

„(°)2 



T+ 



12 



d£e 



f + r (* +21 (%)r + 0[(* - *,-)" 



(k — fcj) 2 y 

Substituting this result into ( |46] ) , we arrive at the perturbation- induced evolution of kj : 



5% _ (o) 



(47) 



We conclude this section with the note that the general (i. e., with account for the perturbation) evolution 
equations for the spectral data 0(r,k), kj(r), and cf\r) are given by addition of their perturbation-induced 



evolution (represented by the variational derivatives (|38|), (|44|), and (|47]) to the evolution equations of these 
quantities in the integrable limit. 
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6 Conclusion 



We have shown that the Riemann-Hilbert problem is a natural setting for dealing with both integrable and nearly 
intcgrable nonlinear equations with the singular dispersion relations. In particular, for small perturbations of 
the integrable models, we have explicit perturbation- induced evolution equations for the RH data, which are 
necessary for analysis and approximate solution for the physically interesting quantities. 
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